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In this paper the geometric entanglement (GE) of systems in one spatial dimension (ID) and in 
the thermodynamic limit is analyzed focusing on two aspects. First, we reexamine the calculation 
of the GE for translation-invariant matrix product states (MPSs) in the limit of infinite system 
size. We obtain a lower bound to the GE which collapses to an equality under certain sufficient 
conditions that are fulfilled by many physical systems, such as those having unbroken space (P) 
or space-time (FT) inversion symmetry. Our analysis justifies the validity of several derivations 
carried out in previous works. Second, we derive scaling laws for the GE per site of infinite-size 
ID systems with correlation length ^ 3> I. In the case of MFSs, we combine this with the theory 
of finite-entanglement scaling, allowing to understand the scaling of the GE per site with the MFS 
bond dimension at conformally invariant quantum critical points. 



I. INTRODUCTION 

Quantum many-body systems in one spatial dimen- 
sion (ID) have proven relevant in physics. For instance, 
quantum spin chains have helped to better understand 
quantum phase transitions [1, 2], the renormalization 
group [ '] , quantum state transfer [4] , and even the fun- 
damental limitations of classical and quantum comput- 
ers [5]. In this respect, the advent of White's Density 
Matrix Renormalization Group (DMRG) [ ] was a major 
breakthrough since it allowed to simulate many of these 
systems efficiently. At its core, DMRG is based on rep- 
resenting physically-meaningful quantum states in terms 
of Matrix Product States (MPSs) [ ]. Thus, the impor- 
tance of MPSs can not be overemphasized: they seem 
to be the correct low-energy description of the relevant 
quantum states of Nature in ID [ ]. 

Also, in recent years there has been growing interest 
in understanding the properties of entanglement in ex- 
tended systems [.)]. While many studies have focused 
on bipartite measures such as the entanglement entropy 
and single-copy entanglement [2], there has also been 
raising interest in investigating the behavior of multi- 
partite quantum correlations. In this respect, among 
many approaches [!)], the so-called geometric entangle- 
ment (GE) [ ] has been demonstrated to be useful 
in a variety of situations, including ID systems and 
MPSs [11-13], 2D systems [ and fully-connected sym- 
metric systems [ ]. However, and as we shall see here, 
the study of GE for ID systems deserves a more detailed 
analysis in some cases. 

The goal of this paper is to analyze in detail two fun- 
damental aspects of the GE of ID systems in the ther- 
modynamic limit. First, we reexamine the calculation of 
the GE of MPSs in the limit of infinite size, offering a 
wide perspective of the problem. As a result, we produce 
a non-trivial lower bound to the GE. This lower bound 
collapses to an equality under certain sufficient condi- 
tions which we shall make precise later. Moreover, we 



show that some of these conditions are natural in many 
physical systems, such as those having unbroken space 
inversion or parity (P) symmetry or the larger spacetime 
inversion (PT) symmetry [. ]. Second, we derive scal- 
ing laws for the GE per site of infinite-size ID systems 
with finite but large correlation length. Also, in the case 
of MPSs we relate these scaling laws with the theory of 
finite-entanglement scaling [ ; . ] and obtain scaling laws 
with the MPS bond dimension at conformally invariant 
quantum critical points. In contrast to the diverging be- 
havior of the entanglement entropy of with the subsytem 
size [ the GE per site is always bounded, and this gives 
certain advantage in terms of convergence in numerics. 

The remaining of the paper is organized as follows. 
In Sec. II we review the definition of the geometric mea- 
sure of entanglement, or in short geometric entanglement 
(GE). In Sec. Ill we discuss the bounds on GE per site in 
the thermodynamic limit. Under certain symmetry, the 
provided upper bound becomes the exact expression of 
GE per site. In Sec. IV we investigate the scaling of GE 
from three different perspectives: (a) finite-correlation 
length, (b) finite bond-dimension, and (c) finite size. We 
summarize in Sec. V. We relegate certain detailed discus- 
sions on the scaling to Appendix A in order to keep the 
smooth flow of the main text. 



II. GEOMETRIC ENTANGLEMENT IN A 
NUTSHELL 

Let us quickly remind the basics of GE. Imagine that 
we are given a quantum state Y^) of TV parties belong- 
ing to a Hilbert space H = {g)f^i VM, where VM is the 
Hilbert space of party r. This could be, for instance, 
the state of a system of quantum spins placed on a lat- 
tice where each party is either a single spin or a block of 
spins. Our aim is to focus on the closest normalized prod- 
uct state of the parties to | ^P) . By "closest" we mean the 
normalized product state 1$) = ((j)^'^'^) ® ■ ■ ■ 
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that minimizes the squared distance 1 1 1$) — j^*) | p between 
1$) and jvl/) or, in other words, maximizes the absolute 
value of their overlap [10], Aniax(*) = max$ K*!^')!. No- 
tice that finding the maximizing product state in this 
equation is actually akin to a "mean-field" approxima- 
tion to I*). 

The larger Amax is, the less entangled is Thus, 
the closest product state approximation to |^') allows us 
to quantify its entanglement via the extensive quantity 
= — log A^a,x(^) where we have taken the nat- 
ural logarithm. As required, E(^) is zero if |^) is unen- 
tangied. We can also define the thermodynamic quantity 
£ and its finite-size version £n as 

£= lim ^AT, £n = N-'^E{^'). (1) 

The quantity £ in the above equation defines the global 
geometric entanglement per site, or density of global geo- 
metric entanglement. This will be the quantity of interest 
in this paper. 

III. BOUNDS OF GE PER SITE IN THE 
THERMODYNAMIC LIMIT 

MPSs offer an accurate description of many interesting 
states of ID quantum many-body systems. For a system 
of size A'^ with periodic boundary conditions, these are 
states defined as 

1*)^ ^ tr(AW^^...A[^l^~)|zi,...,i^), (2) 

ii,...,ijv 

where AI™!*'" is a x x X niatrix at site m — 1, . . . ,N for 
each i„i = 1, . . . which labels a local basis of the 
Hilbert space of dimension dm at site m. Parameter x is 
called the bond dimension of the MPS. It is known that 
any spin state can be written as a MPS, albeit the bond 
dimension may need to increase with the system size [ ]. 
Therefore, MPS are a good variational family of states to 
approximate any state, and in particular, ground states 
of Hamiltonians with local interactions. If the system is 
invariant under translations of one site (TI) [_ >], then 
we can assume to have the same matrices at every site, 
j\[m]i^ _ This allows us to take the thermo- 

dynamic (infinite-size) limit of an MPS with TI just by 
considering the same matrices A^"^ at all the infinitely- 
many sites m = 1, . . . , oo (see e.g. Refs. [20]). 

Let us now discuss the calculation of the GE of an MPS 
in the thermodynamic limit. One important assumption 
that we make now is that the closest product state |<I>) 
can be taken to be a product of identical local states, that 
is, it fulfills TI: |$) = |<^)«'°°. This happens to be a good 
choice for e.g. the ground states of the transverse-field 
XY spin chain, where this ansatz has been verified numer- 
ically [11]. For permutation invariant states, this ansatz 
has been proven correct [17]. However, we caution that 
sometimes this assumption does not hold. For instance, 
the state l^-) = 2-i/2(|0101...) -|- [1010.. .)) is manifestly 



TI, whereas its closest product states (e.g. [0101. ..) or 
1 1010...)) are not. In such cases a modified product state 
ansatz with alternating periodicity (e.g. TI every 2 sites) 
is the proper choice, as has been verified numerically [ 1 2- 
14, 21]. Moreover, one can instead consider GE defined 
with respect to block product states, and as long as the 
size of a block is chosen to be a multiple of possible peri- 
ods of the possible symmetry-breaking states, such trans- 
lation invariant ansatz is expected to hold. 

Under the above assumptions, the overlap between an 
infinite MPS \'^) and a product state |$) reads 

($[*) = hm tr(i?^), (3) 

where is a x x x transfer matrix defined as = 
j4*((/)|j). Let us call the eigenvalues of matrix B^ 
(where a labels the different eigenvalues). Then, we have 
that tr(B^) = J2J\'})^ - fc(Aj,)^, where the last ap- 
proximation is valid in the limit N ^ oo and where we 
assume Aj, to be the eigenvalue of largest absolute value 
and k its total degeneracy. 

The absolute value of the eigenvalue Aj, is also known 
as the spectral radius of matrix B^, = [A^]. Another 
important quantity is the numerical radius of i?^, defined 
as Wcj, = maxf I B^ r |, where || r || = 1. A crucial 
property is that the numerical radius of a matrix upper 
bounds its spectral radius [23], 

P<f, <w^. (4) 

The above expression turns into an equality if (but not 
only if) matrix B^ is diagonalizable. 

Our bound for the GE is based on considering in- 
stead of p0 as the relevant dominant scale in the thermo- 
dynamic limit for the overlap in Eq. (3). As the latter is 
not larger than the former, we have specifically that the 
maximum overlap over product states is bounded as 

Amax(*) ^ max{p^) < max{w^) ( max(w0) ) (5) 
4> <P \ <t> J 

where again the first approximation is valid when N 
tends to infinity (and up to an irrelevant degeneracy fac- 
tor independent of A^). Therefore, in the thermodynamic 
limit the density of geometric entanglement £ in Eq. (1) 
obeys the bound £ > — 2 log (max^ maxf | B^ ^\)- 
The advantage of the latter expression is that it is vari- 
ational. Furthermore, the double maximization in the 
right hand side of this equation can be further simpli- 
fied under the assumption that both maximizations com- 
mute [12, l"]. In such a case, the maximization over |0) 
can be done straightforwardly and we obtain 

£ > -2 log ^max I (f®r*)'f E (r«)f *)|^ , (6) 

where E is the zero-dimensional ^ MPS transfer 
matrix E = ® {A^)* . The above equation is our 

main lower bound on the density of GE. 
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Let us now discuss the result in Eq. (6). To start with, 
notice that if the MPS is such that is always di- 
agonalizable V |(/)), then the spectral radius and the 
numerical radius coincide and Eq. (6) turns into an 
equality for the density of GE. Notice also that this is 
a sufficient condition, but not necessary. In fact, the 
"equality" version of Eq. (6) was originally derived in 
the second paper of Ref. [12], where certain assumptions 
guaranteed the diagonalizability of i?^ for any \(j>). More 
specifically, it was the (unbroken) symmetry of the phys- 
ical system under space inversion (P symmetry) the key 
property that guaranteed this diagonalizability. In order 
to see this, consider a translation invariant MPS state 



J2 tr(A*i •••A'")|zi 



(7) 



which is already in Jordan normal form and such that 
its eigenvalues satisfy A J, > A^. For simplicity let us 
assume that these eigenvalues are real and positive. It 
is easy to check that the spectral radius of this matrix 
is given by = A J,, whereas the numerical radius is 
= max(Aj,,A^ + 1/2). According to this, we have 
that for A^ > A^ + 1/2 the numerical and spectral radius 
coincide and then Eq. (6) collapses to an equality if such 
a behavior is found for any state whereas for A^ < 
A| + 1/2 both radius are different and therefore Eq. (6) 
must remain as a lower bound. 



IV. SCALINGS OF GE PER SITE 



and its mirror image state (by reversing the ordering in 
a iV-spin basis state without changing the amplitude) 



ii,...,ijv 

In order for l^f) = \'^'), it is sufficient to have {A")'^ — 
UA'U-\ for any invertible U [24, 25]. Thus, {A')'^ = A" 
is a sufficient condition for the space inversion or parity 
(P) symmetry. Such a condition guarantees that any i?^ 
is automatically symmetric, and therefore always diago- 
nalizable by some orthogonal transformation. P (or the 
larger PT) symmetry of the physical system is thus a suf- 
ficient condition for any matrix to be diagonalizable, 
which in turn is a sufficient condition for Eq. (6) to col- 
lapse to an equality. We stress, though, that a matrix 
B^ may still be diagonalizable even if the MPS does not 
have P (nor PT) symmetry. 

The general case may be more intricate. To see this, 
notice that any matrix B^ (diagonalizable or not) can 
be written as = P^^J^P^, where is the Jordan 
normal form of B^ [-■ ]■ Matrix is a direct sum of 
Jordan blocks, namely = 0^ -^^ with {J^)i^,n = 
A^^i,™ + A^ being the eigenvalues of B^, and 

l,m = 1, . . . , go with qa the size of the Jordan block a. 
An important theorem in linear algebra states that a ma- 
trix is diagonalizable if and only if gc, = 1 Va, that is, all 
Jordan blocks are trivial. Thus, the appearance of non- 
trivial Jordan blocks in the Jordan normal form of B^ 
makes it non-diagonalizable, which in turn implies that 
Eq. (6) may need to remain as a lower bound and not 
collapse to an equality, depending on the properties of 
the Jordan blocks in J^. In such cases, the true value 
of the GE must be computed using the left hand side of 
Eq. (5) and maximizing the spectral radius over [</)), 
which may be non-trivial. For instance, consider the case 
of the non-diagonalizable matrix 



B, 




(8) 



In this section we investigate the scaling of GE from 
three different perspectives: (a) finite-correlation length, 
(b) finite bond-dimension, and (c) finite size. 



A. Finite-^ scaling of GE per site 

The GE also obeys precise scaling relations (see e.g. 
Ref. [12]). We now focus on analyzing the scaling of the 
density of GE in Eq. (1) as a function of the correlation 
length ^ for infinite-size ID systems, where we assume the 
correlation length to be large. Generally speaking, since 
the GE per site is a density quantity (similar to the free 
energy density) the singular part of GE near criticality 
is expected to behave as 

f (0 - £(e = oo) ^ (6 -(- &' log^ + • • • (9) 

where d is the spatial dimension, and 5, h' model- 
dependent coefficients. In this expression, we have in- 
cluded a first non-trivial logarithmic correction and as- 
sumed that the omitted parts are less singular. In order 
to assess the validity of Eq. (9) we consider some ana- 
lytical examples, namely (i) the case of MPSs for which 
Eq. (6) turns into an equality, and (ii) the quantum XY 
model for spin 1/2 in the Ising and XX regimes. 

Let us then start with the case of MPSs. Under certain 
assumptions, it is possible to find the corrections due to 
finite correlation length for the GE per site of an MPS for 
which the equality version of Eq. (6) holds. First, let us 
assume that the MPS transfer matrix E obeys an spec- 

tral decomposition E = X]q=i A*"^"-^"^; where i?" and 
L" are respectively the ath right and left eigenvectors of 
E with eigenvalue /x". If the eigenvalues /i" are different 
and rapidly decaying, then this spectral decomposition 
can be approximated hy E ^ ^^{R^L^^ + {^'^ / ij})B?L'^^), 
and therefore £ ^ — 21oglp -I- {p? / ij})q\, with p and 
q some given coefficients. At this point we recall the 
standard definition of the correlation length of an MPS, 
^ = — 1/log l^^/^^J. Using this, together with the as- 
sumption ^ ^ 1, it is not difficult to arrive at the ex- 
pression £{^) ^ a + b/£^, with a — — 21og]p -|- q\ and 
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b = 2q/\p + q\. We see thus that an MPS is in principle 
capable of handling 0(1/^) corrections. Nevertheless, al- 
ternative scaling relations may also hold for MPSs if some 
of the considered assumptions break down. For instance, 
it could be the case that the eigenvalues of the MPS 
transfer matrix E are not rapidly decaying, and there- 

2 

fore £ ^ ~2 log I ^^^]^ /i"(7"| for the GE (where q" are 
some coefficients). In principle, this more general expres- 
sion could account for all the terms in the general scaling 
law from Eq. (9). 

Next, we consider the case of the ID XY model for 
spin 1/2, defined by the Hamiltonian 



(a) 



H 



XY 



E 



1 



1 



hal: 



(10) 

where h is the magnetic field, r is the anisotropy pa- 
rameter, and (Tq' is the ath Pauli matrix at site i. For 
< r < 1 the model belongs to the universality class of 
the ID quantum Ising model (central charge c = 1/2). 
In this regime, in Ref. [ ] it was proven that 



dh 



1 



2'Kr log 2 



log|/i- 1|, for 1| < 1. (11) 



Using the above equation together with the known result 
^ = l/\h — 1|, we can obtain the dependence of £ on ^. 
Integrating the expression in Eq. (11) we obtain 



(12) 



for ^ ^ 1 and where a, b and b' are some coefficients 
independent of Notice that this is the same type of 
scaling as in Eq. (9). Also, for r = the model belongs 
to the universality class of the ID quantum XX model 
(central charge c = 1). In this case, and near h = 1^ , 
the GE per site was shown [11] to behave as £{£,,r — 
0) 1/C, where ^ ~ |1 - h\-^^^. It should be pointed 
out that although the XX model near this transition is 
scale invariant, it is not conformally invariant. However, 
near h — we have seen numerically that a law like 
the one in Eq. (12) seems to emerge also for the XX 
model (see Appendix). Remarkably, we see that all these 
scaling behaviors accommodate as well in the law hinted 
in Eq. (9). 
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FIG. 1: Exponent Ki cLS Si function of (a) r for the XY model 
with h — 1, and (b) h for the XX model. The values and 
errorbars correspond to those in tables I and II. Dotted lines 
correspond to the averages over the studied interval, k ~ 2.1 
for (a) and k ~ 1.3 for (b). 



quantum critical point the corre lation length obeys the 
relation ^ ~ x'^ with k = 6/ {c{y/l2/c + 1)) (c being the 
central charge of the associated universality class) [19], 
one obtains the leading term expression 



£{x}-^a+{b + b'\ogx)/x^ 



(13) 



with a, b and b' coefficients. Remarkably, the validity of 
this relation is in agreement with numerical simulations 
with MPSs for a variety of ID spin chains, within some 
accuracy considerations (see Appendix). As illustrated 
in Fig. 1, the extracted exponent n is very close to the 
expected value. We suspect that the deviation is due to 
the fact that (i) the GE converges quickly with x to the 
X = oo value and (ii) the correction is small (as it scales 
inversely with x to some positive power). The larger 
deviation for smaller r in Fig. la and that at larger h in 
Fig. lb is also due to the fact that the r — and h — 1 is 
not a conformally invariant critical point, and the scaling 
does not hold. For further discussions, see the Appendix. 



C. Finite-size scaling of GE per site 



B. Finite-x scaling of GE per site 

An important fact is that the finite-^ scalings con- 
sidered here can be combined with the theory of finite- 
entanglement scaling [ ] , in turn allowing to understand 
the finite-x scaling of the GE per site of an MPS with 
bond dimension x- In contrast to the diverging behavior 
of the entanglement entropy of with the subsytem size [2] , 
the GE per site is always bounded above by log2 d ebits, 
where d is the local Hilbert space dimension, and this 
gives certain advantage in terms of convergence in nu- 
merics. Knowing that close to a conformally-invariant 



Finally, our scaling laws for finite correlation length are 
in accordance with the finite-size scaling behavior found 
for ID systems of size N at the thermodynamic critical 
point [12, 13, 21, 22]. In that case, it was found that 
the GE per site obeyed the law £{N) ^ a + b/N, with 
a and b some size-independent coefficients. Up to a log- 
arithmic factor in some cases, this is the same type of 
behavior that we have found here if the role of the cor- 
relation length ^ is now played by the size of the system 
N. In general, this is a manifestation of the well-known 
property that for a finite system close to criticality, the 
size of the system plays the role of an infrared cut-off in 
the correlation length. 
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V. CONCLUSIONS 

Here we have provided a non-trivial lower bound to the 
GE of MPSs in the thermodynamic limit. We have dis- 
cussed some sufficient conditions under which this lower 
bound collapses to an equality, such as space inversion or 
parity (P) symmetry, in turn justifying the approach con- 
sidered in previous works to compute the GE. We have 
also derived scaling laws for the GE per site of infinite- 
size ID systems with finite but large correlation length 
^ 1. These results have also been related to scaling 
of the GE with the bond dimension of MPSs and the 
finite-size scaling of the GE for ID systems. 
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Appendix A: Further discussions on finite-x scaling 

Here we discuss numerical results that further check 
the validity of Eq. (13). We have performed numerical 
simulations for the spin-1/2 quantum XY model in the 
regimes corresponding to both the Ising and XX univer- 
sality classes (respectively central charges c = 1/2 and 
c = 1), and for the spin-1/2 antiferromagnetic Heisenberg 
model (c = 1). The MPS approximation to the critical 
ground state has been obtained using the iTEBD method 
[20] for X G [2,40], and the GE has been extracted using 
standard optimization tools. 

Before presenting our results, a word is in order re- 
garding potential sources of errors. As opposed to other 
quantities such as the entanglement entropy or the single- 
copy entanglement, the GE at critical! ty converges fast 
to a finite value as x (the MPS bond dimension) grows. 
This is indeed a good feature if we are interested in the 
GE itself in the limit of infinite bond dimension (as op- 
posed to e.g. the entropy, which tends to be divergent). 
Our fits try to capture the correction to this converged 
finite value. In the vast majority of the cases, such a 
correction is small and approaches zero as x — >■ cx) very 
quickly. Thus, our fits are sensitive to small numerical 
errors difficult to control completely. These errors are ac- 
centuated if, on top, the GE itself is very small, as is the 
case of the XY model close to r = and h = 1 (where 
the ground state of the system is separable). If this is 
the case, then our fits try to capture a tiny correction 
to a tiny quantity, which may be numerically ill-defined. 
As we shall see, this seems to be particularly true for the 
XY model at ft. = 1 and as a function of the anisotropy 
r, where the GE itself is around one order of magnitude 
smaller than for the XX and Heisenberg models. In turn, 
this also implies that the our most accurate fit to Eq. (13) 



r 


a 




b' 


b 


0.5000 


0.015 


2.76 ±0.11 


0.010 ±0.001 


-0.017 ±0.002 


0.5625 


0.017 


2.50 ±0.14 


0.012 ±0.001 


-0.020 ±0.002 


0.6250 


0.019 


2.37 ±0.15 


0.013 ±0.001 


-0.023 ±0.002 


0.6875 


0.021 


2.18 ±0.14 


0.014 ±0.001 


-0.026 ±0.002 


0.7500 


0.023 


2.07 ±0.15 


0.014 ±0.001 


-0.028 ±0.003 


0.8125 


0.024 


2.00 ±0.15 


0.015 ±0.002 


-0.030 ±0.003 


0.8750 


0.026 


1.93 ±0.15 


0.016 ±0.002 


-0.032 ±0.003 


0.9375 


0.028 


1.89 ±0.15 


0.017 ±0.002 


-0.034 ±0.004 


1.0000 


0.029 


1.85 ±0.15 


0.017 ±0.002 


-0.035 ±0.004 



TABLE L Value of the fitted parameters a, b, b' and the 
exponent k together with the associated fitting errors, for the 
XY model with h — 1, for different values of r. The fits are to 
Eq. (13), and in the regime where the GE is not exceedingly 
small. On average, k ~ 2.1. Parameter a has no fitting error 
since it is extracted directly from our numerical algorithm in 
the X ~^ oo limit. 



is obtained for the Heiseiberg model, for which the GE 
is larger. Nevertheless, and in spite of these considera- 
tions, our fits succeed in capturing the essential scaling 
properties of the different systems with good confidence 
in some regimes. 

Keeping the above considerations in mind, our numer- 
ical analysis indicates the following: 

(i) For the XY model with r 7^ and h — 1 (central 
charge c = 1/2 universality class), the law in Eq. (13) fits 
well the GE for r close to 1, whereas for r far from 1 the 
GE decreases quickly towards very small values. This 
induces to variations in the different parameters of the 
fit as a function of r, which become stronger as the GE 
approaches zero (see the data in Table I). The observed 
variations in the parameters are compatible with the laws 
a{r) ^ aa + Par, b{r) -ab - Pbr, b'{r) ^ ab' - /r^/^ 
and K{r) ~ ± /J^/r^, for some positive coefficients 
Ota, PajOtb, j3b,0Lb' Pb' ,OLi^ and /3„. In spite of this varia- 
tions, the average behavior of the scaling exponent seems 
to be K '--^ 2.1, compatible with the theoretical prediction 
K - 2. 

(ii) For the XY model with r = and < ft < 1 
(central charge c = 1 universality class) , we obtain good 
fits to Eq. (13) in the region close to ft = with values 
of the scaling exponent around k ^ 1.2 (see the data in 
Table II). Close to ft = 1, where the ground state be- 
comes very little entangled, the GE is again too small 
to extract the scaling behavior with confidence. Further- 
more, the critical point near ft = 1 is not conformally 
invariant and has dynamical exponent z = 2, and hence 
the relation between k and c (the latter being not de- 
fined in this regime) does not hold. The observed vari- 
ations in the parameters are compatible with the laws 

a(ft) ^ aa- l3ah^, h{h) ab ± /36ft^ b'{h) ~ ab' - Pb'h 

and niji) ^ ± /3Kexp(7Kft), for some positive coeffi- 
cients aa,l3a,ab,l3b,ab'Pb',a^,l3i^ and 7^. Interestingly, 
we obtain /3k '-^ 3 x 10~^ and 7^ ^ 10, which involves 
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h 


a 




b' 


b 


0.0000 


0.155 


1 


21 ±0.06 


0.060 ±0.003 


-0 


171 ±0.008 


0.0625 


0.154 


1 


21 ±0.06 


0.060 ±0.003 


-0 


171 ±0.008 


0.1250 


0.154 


1 


21 ±0.06 


0.060 ±0.003 


-0 


170 ± 0.008 


0.1875 


0.152 


1 


21 ±0.06 


0.059 ±0.003 


-0 


168 ± 0.008 


0.2500 


0.151 


1 


21 ±0.06 


0.059 ±0.003 


-0 


167 ±0.008 


0.3125 


0.148 


1 


22 ±0.06 


0.059 ±0.003 


-0 


166 ± 0.008 


0.3750 


0.145 


1 


22 ± 0.06 


0.057 ±0.003 


-0 


162 ± 0.008 


0.4375 


0.141 


1 


22 ±0.07 


0.057 ±0.003 


-0 


159 ±0.008 


0.5000 


0.137 


1 


23 ±0.07 


0.056 ±0.003 


-0 


156 ±0.009 


0.5625 


0.132 


1 


23 ±0.07 


0.055 ±0.003 


-0 


152 ±0.009 


0.6250 


0.125 


1 


25 ±0.08 


0.054 ±0.004 


-0 


146 ± 0.009 


0.6875 


0.118 


1 


29 ± 0.09 


0.054 ±0.004 


-0 


142 ±0.010 


0.7500 


0.108 


1 


33 ±0.10 


0.053 ±0.005 


-0 


136 ±0.011 


0.8125 


0.097 


1 


40 ±0.12 


0.051 ±0.006 


-0 


128 ±0.012 


0.8750 


0.081 


1 


64 ±0.21 


0.061 ±0.011 


-0 


133 ±0.021 



TABLE II: Value of the fitted parameters a, b, b' and the 
exponent k together with the associated fitting errors, for the 
XX model for different values of h. The fits are to Eq. (13), 
and in the regime where the GE is not exceedingly small. 
On average, k ~ 1.3. Parameter a has no fitting error since 
it is extracted directly from our numerical algorithm in the 
X — >■ 00 limit. 



a smooth behavior of the exponent for low values of h, 
with a rapid increase for large h. Despite of this varia- 
tions, the average scaling exponent that we obtain in this 
region is k ~ 1.3, compatible with the theoretical result 
(which is also k ~ 1.3). 



(Hi) For the spin- 1/2 Heisenberg model (central charge 
c — I universality class), our fit to Eq. (13) seems to indi- 
cate that a - 0.259, b 0.259±0.004, b' - 0.078±0.010 

and K ^ 1.25 ±0.08, compatible also with the theoretical 
result K ~ 1.3. 



Thus, our data seems compatible with the scaling law 
in Eq. (13) with values of the scaling exponent n not 
too far from the theoretical predictions for conformally- 
invariant quantum critical points, within the considera- 
tions of numerical accuracy that we discussed previously. 
Let us also remark that, independently of the specific 
form of scaling, the fact that the GE converges quickly 
with X to a specific value is also a remarkable and use- 
ful property, as e.g. the infinite-x limit can be easily 
extracted. 
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